Testing quantum correlations in a confined atomic cloud by scattering fast atoms 
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We suggest measuring one-particle density matrix of a trapped ultracold atomic cloud by scatter- 
ing fast atoms in a pure momentum state off the cloud. The lowest-order probability of the inelastic 
0^ , process, resulting in a pair of outcoming fast atoms for each incoming one, turns out to be given 

ON . by a Fourier transform of the density matrix. Accordingly, important information about quantum 

correlations can be deduced directly from the differential scattering cross-section. A possible design 
O ■ of the atomic detector is also discussed. 
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Successful advances in achieving collective quantum states in confined clouds of alkaline atoms pj| and in atomic 
hydrogen [^j make possible studying quantum coherent properties of these systems, as well as the revealing fundamental 
kinetic processes leading to the formation of the coherence. In Ref. ||, a dynamics of the condensate growth has been 
observed by means of detecting a density increase usually associated with the condensate formation in the trap. The 
coherence of the condensate has been demonstrated directly by letting two released condensates form the interference 
pattern f|. 

Mechanism of formation of quantum correlations is a matter of great attention and controversy. An emergence 
' of the condensate due to either lowering temperature or reaching an equilibrium after fast quenching is associated 
with a formation of the off-diagonal long-range order (ODLRO) §. A primary object displaying such an order is the 
one-particle density matrix (OPDM) p{x.\, x 2 ). Typical distances over which these correlations become important are 
O , comparable with the interatomic spacing r a . Consequently, an "early" detection of such emerging correlations is very 
difficult to achieve by light with the wavelength A>r a usually employed for probing the cloud density. In Ref. ||, 
t-H , it has been suggested that a resonant fluorescence of an external atom identical to those forming the condensate is 
sensitive to the relative phase of two condensates. However, this method is not suitable for testing distances shorter 
than A. Information about short-range density correlations (at distances < r a ) can, in principle, be obtained from the 
absorption of detuned resonant light 0. The change in local m-body density correlations (the so-called m!-effect Q]) 
, can be seen by measuring recombination rates, and this has been already done experimentally for the equilibrium case 
|)| . However, measurements of the correlation length r c of the forming ODLRO seem very unlikely to be achieved by 
these methods. 

Thus it is tempting to have a tool which could make possible seeing the OPDM directly without limitations on the 
accessible distances. 

Currently, great efforts are being dedicated to a creation of a controllable source of coherent atoms - atomic LASER 
(see Ref. |l(} and references therein). In Ref. jDJ, a mechanism of accelerating neutral atoms has been proposed. 
Hence, it is very likely that a source of fast and coherent atoms will be available soon. In this paper, we suggest a 
method of detecting the OPDM which relies on inelastic scattering of such atomic beam off the atomic cloud. 

We note that methods of scattering of neutrons |1^|l|] and He atoms |l4|] off liquid He are well known. In Ref. 
lyl, the Impulse Approximation has been suggested to employ for interpreting the differential cross-section of fast 
neutrons. In such an approximation, it is possible to relate the momentum transfer distribution to some integral of 
the population factor in He. However, liquid He is a strongly interacting system where the gas parameter £ = na 3 
(n is the density and a is the scattering length) is not small. Accordingly, multiparticle excitations dominate in the 
final-state channel, which makes a direct measurement of the OPDM impossible. Certain assumptions about the role 
of the final-state effects should be made Jl3|] . In contrast to this, the gas parameter in the trapped atomic condensates 
can be as small as £ ~ 10~ 5 . Accordingly, the mean free path is « a/£ ~ 10~ 2 cm, which greatly exceeds the mean 
particle separation, and can become greater then the cloud size. This implies that a contribution of the multiparticle 
excitations can be safely neglected as long as a wavelength of the incoming atom is much smaller than the interparticle 
separation in the atomic cloud. Under these conditions, as we will discuss below, it becomes possible to measure the 
OPDM directly. 

If the mass of the incoming fast atom is comparable to that of the atoms forming a cloud (or, in particular, the 
fast atom is just identical to the atoms of the cloud), the lowest-order inelastic scattering event is a production of two 
fast outcoming atoms (the size of atomic cloud is supposed to be small enough to neglect multiple scattering). The 
quantum-mechanical probability of this process turns out to be proportional to a Fourier transform of the OPDM. 
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Another process - an elastic scattering of a single fast atom off the cloud - also exists. In this case one has a single fast 
atom in the final state, in a complete analogy with elastic scattering of light. In principle, the elastic processes might 
mask the inelastic ones. However, if the momentum of incoming atom, k, is much larger than the typical momentum 
of the the cloud, ~ l/r c , then the elastic-scattering angles are small in the parameter (fcr c ) _1 <C 1. This feature 
allows one to distinguish between the two processes. One more feature, arising at kr c ^> 1 from the conservation laws, 
is the fact that the angle between the two created fast atoms is always close to ir/2. This facilitates identification of 
pairs resulted from one and the same scattering event. 

Apart from extremely small uncertainty of the momentum of the incident atom, the method we suggest implies 
that the total momentum of the outcoming pair can be detected with a sufficient precision. Thus, a reliable atomic 
detector is required. Therefore, we will also discuss a possible design of such a detector. 

Now let us derive an expression for the inelastic cross-section in the lowest order with respect to the two-body 
interaction. For the sake of simplicity, we will not consider the complexity of all possible scattering channels and will 
concentrate on a simplest case of spin-polarized bosons, when the incoming fast atom is identical to particles forming 
the cloud, and its spin polarization is the same as that of the cloud. We emphasize that the validity of the suggested 
method relies on a possibility to have the incident momentum k obeying the relation 

£ 1/3 <fca<l. (1) 

Then, one can represent the interaction Hamiltonian in a traditional form (h = 1) 

H int = ^ / dx*t( x )tft( x )tf( x )tf( x ) UQ = ^ (2 ) 
2 J to 

with to standing for the atomic mass. The total field 'J can be subdivided into the low- and the high-energy parts, 
ip and ip' , respectively: 

* = ^ + V', V' = 5> k e lkx , (3) 

k 

where the incident and the scattering states are described in terms of plane waves normalized to unit volume; at 
destroys the high-energy particle with the momentum k. A substitution of Eq. (||) into Eq. (|J) and selection of the 
terms that describe a process involving one incident fast atom carrying momentum k and two ejected fast atoms 
carrying momenta ki and k2, as well as its reverse, yield 

H- nt = u J2 a kl a k 2 fl k / rfxe " q ^ + H.c. , (4) 

k!,k 2 ,k 

where q = ki + k2 — k is the transferred momentum. The double-differential cross-section for scattering with given 
energy and momentum transfers, to and q, in the lowest order with respect to H' int is given by the Golden rule formula 
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W(q,«) = J ^ &[u,-u> fi ) J J d Xl d X2 J die^ xl - X2 )--V(xi,t;x 2 ,0) , 



(5) 



where /o(xi, t\\ x 2 , £ 2 ) — (^(xi, ii)^^, f 2 )) [Note that unlike the uniformity in time, p{x-i, t\; x 2 , i 2 ) = /°( x i^i — 
i 2 ;x 2 ,0), the space uniformity cannot be, in general, assumed as long as a trapping potential exists.]; ojfi is the 
difference of the kinetic energies of the fast atoms in the final and initial states: 



UJfi = — 
TO 
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— + qk + k\ - (q + k)k x 



(6) 



with ki being the momentum of one of the outcoming atoms while the momentum of the second one, k 2 , is fixed by 
the relation k 2 = q + k ki. 

The integration over ki in Eq. (^) can be carried out explicitly. However, first we notice that the requirement of 
large k means that the values of q and u>, which are effectively selected by the correlator p(xi,t;x 2 ,0) in the right- 
hand side of Eq. (ph, satisfy the conditions q <C k and \u>\ <C k 2 /m. This immediately leads to the approximation 
5(u> — uifi) Ri TO(5(fcj — kki), which yields 

W{q,tu) =4a 2 J J dxidx 2 J dte^ 1 '^-^ p(x lf t; x 2 , 0) . (7) 
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We thus see that the double-differential cross-section W(q, u>) is almost directly related to the dynamic correlator 
p(xi, i; x 2 , 0) which, as is known, contains rather rich information about the system, including, for one thing, the 
elementary excitation spectrum. Confining ourselves to the static correlations, described by OPDM / o(x 1 ,X2) = 
p(xi, 0; X2, 0), we arrive at even more simple relation 

MK(q) = 87ra 2 / / dxidx 2 e ,;q < Xl - X2 V(xi, x 2 ) (8) 



in terms of the differential cross-section W(q) — J duj W(q, uo). 

Consider the structure of VF(q) in the most characteristic cases. In the case of pure Bose-Einstein condensate, 
when p(xi,x 2 ) = $*(xi)$(x 2 ) [^(x) is the condensate wave-function], 

W(q) = 87ra 2 |$ q | 2 , $ q = J dxe-' iqx $(x) . (9) 

In particular, if <f> q is real, Eqs. @ can be reversed, and correspondingly $(x) can be restored from M^(q). 

Another instructive case is the axially symmetric quantum vortex. A presence of a single vortex in a center of the ax- 
ially symmetric condensate drastically changes the scattering pattern. Indeed, in this situation, $ = cnp(i8)y / 'n(r, z), 
where 9 is the axial angle and n(r, z) stands for the axially symmetric condensate density as a function of the dis- 
tances r, z perpendicular to the axis and along the axis, respectively. Accordingly, Eq. (||) shows that ^(q) = for 
q directed along the vortex line. The differential cross-section becomes finite as long as there is a component qj_ of q 
perpendicular to the axis or the vortex displaces from the condensate center. In the first case, W(q) ~ q\ for q± — > 0. 

Quite similar to the quantum vortex is the case of the supercurrent state of a toroidal Bose condensate. The 
suppression of W(q) for q perpendicular to the plane of the torus allows distinguishing the supercurrent state from 
the currentless genuine ground state. 

In the absence of condensate, or for the above-the-condensate part of OPDM, one normally deals with the so-called 
quasi-homogeneous regime, when typical inverse momentum of particles is much less then the scale of density variation. 
In this case, it is reasonable to introduce the variables r = x 2 — X! and R = (xi + x 2 )/2, and to represent W^(q) as 

W(q) = 8na 2 J dRp q (R) , p q (R) = / dxe _iqx p(x, R) , (10) 

since in the quasi- homogeneous regime, OPDM in the Wigner representation, p q (R), has a semiclassical meaning of 
local (at the point R) distribution of the particle momentum q. We thus see that in this case W^q) yields spatially 
averaged momentum distribution. Moreover, without contradiction with the uncertainty principle, this averaging can 
be partially (totally) removed by collimating the incident (both the incident and outcoming) beams. Even without 
removing the averaging, VF(q) contains valuable information about long-range correlations in the system, since the 
the averaging does not affect the order-of-magnitude value of the correlation radius r c (equal to the inverse typical 
momentum k c ). If the system contains both condensate and quasi-homogeneous above-the-condensate fraction, the 
cross-section W^(q) is given by a combination of Eqs. (^) and (|To)) . 

Now let us discuss one possible method of detecting a total momentum of an outcoming pair of atoms. There is a 
significant requirement for such a detector: an uncertainty of detecting the momentum transfer must be less than the 
inverse correlation length r c in the cloud. This imposes a limitation on the detector size Rjj. Indeed, let us suppose 
that cells sensitive to the arrival of the scattered atoms are located on a sphere of the radius Rrj, and the atomic 
cloud of the size R <C Rd is at the center of this sphere. Then, the uncertainty in the scattering angle is R/Rd- This 
produces the uncertainty in the momentum ~ kR/Ro- This uncertainty must be much less than l/r c , or 

R D > kr c R. (11) 

It can be shown that the absolute values of the scattered momenta measured by the time of flight are also subjected 
to the same uncertainty. Thus, the condition ( |lT|) determines a precision of the time-of- flight measurements as well. 

Our method implies that two outcoming atoms produced by one incoming atom can be identified. As mentioned 
above, this identification can be done due to the w 7r/2 angle between the scattered atoms. However, if the incoming 
coherent beam produces too many scattered pairs, the erroneous identification is very likely. Let us derive a limitation 
on the number N s of the scattered pairs per typical time of flight, which would keep the momentum uncertainty less 
than l/r c . The scattering events under consideration should occur with approximately equal probability at any angle 
inside the Air solid angle, so that the angular area occupied by a single event is 4ir/N s . It should be much less than 
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2it/r c k which determines the angular area of the strip where the second atom of a pair can be found, the direction of 
the first being fixed. Thus, 

N s < 2kr c . (12) 

If a typical speed of the fast atoms is lOcm/s, and r c w 10~ 3 cm, R w 10~ 2 cm, the conditions ( [Tl| ) and ( |i~2| ) yield 
Rd ^> 10 cm and N s <C 10 3 . Accordingly, a typical time of flight will be longer than 1 s. 

Now let us discuss a mechanism of detection of a single neutral atom. We suggest employing resonant atomic 
fluorescence in the evanescent field of light propagating inside a waveguide. Specifically, light sensitive cells are 
mounted on a side of a long waveguide so that light propagating inside the waveguide does not excite these cells. If 
holes are made through the cells and the waveguide, the light will remain confined as long as a diameter of the hole 
is much smaller than the wavelength of light. However, a single atom may enter the hole and feel the resonant field 
inside the hole, or in the close vicinity of the entrance to it, due to the evanescent field of the light. Consequently, the 
atom will reemit resonantly one or several photons. These photons can then be detected by the nearest cell indicating 
an arrival of an atom at a specific location at the detector surface. Accordingly, a scattering angle of the coming atom 
can be identified. The velocity of the atom can be deduced from the time of flight. 

A probability Pd for atom to penetrate into a hole is simply the area Ah of the holes per unit area of the detector 
exposed to the atomic flux. In order to detect an atom in the hole, the atom must reemit at least few photon which are 
captured by the cell mounted outside the waveguide. Before we estimate a number Nu of the reemitted photons which 
can be detected, we note that the photons remitted by the atom, which has penetrated deeply inside the waveguide, 
should remain confined inside it, and therefore they will not be detected by the cell. Only those photons, which are 
reemitted by the coming atom while being close to the hole entrance, will be scattered almost isotropically and can 
be absorbed by the cell. Taking into account that the penetration length of the evanescent field is comparable with 
the wavelength A of light, one can find the time to ~ X/v during which the atom moving at speed v is subjected 
to the evanescent field and reemits light isotropically. Then, the number of photons reemitted during this time is 
Nd ~ 7^d = jX/v, where 7 stands for the natural width of the line. Choosing typical values A = 700nm, 7 ps 10 7 
Hz and v = 10 cm/s one finds N p h — 70. Some geometrical factor of the order of one should reduce the number of 
photons reaching the cell. Quantum efficiency of photodetectors can be easily achieved to be 10%-20%. This implies 
that once an atom entered a hole, it will be detected with high certainty. Therefore, a probability to detect an atom 
at a given position is just Pp. In other words, if a diameter of each hole is 100 nm and two closest holes are 300 
nm apart, the probability is Pd ~ 0.1. Thus, in order to achieve a resolution of at least .1 of the typical momentum 
region in the scattering intensity, one needs at least 10 scattering events for each component of the momentum. This 
corresponds to 10 atoms ejected from the cloud (plus 10 fast incident atoms). Accordingly, for the 3D geometry, 
it translates into 10 3 /Pd « 10 4 atoms in total. In a typical condensate of 10 6 atoms such a bombardment by 
fast atoms will result in a depletion of the cloud by only 1%. Evidently, this method is not appropriate for clouds 
containing less than 10 5 atoms. In the effective 2D or ID geometries, the same resolution requires much less numbers 
of the scattering events. Specifically, 10 3 and 10 2 of them for the 2D and ID geometries, respectively, will satisfy the 
above criteria of resolution. 

In conclusion, we have suggested a method of scattering of fast atoms in a pure enough momentum state off a 
trapped atomic cloud in order to test directly one-particle density matrix of this cloud. The differential cross-section 
of the inelastic process, when one incoming fast atom produces two fast ones, allows measuring the correlation length 
of the local off-diagonal order. This gives, in particular, a powerful tool for testing different scenarios of formation of 
the off-diagonal long-range order in the traps. This method can also be employed for detecting quantum vortices and 
supercurrent states, as well as the effect of quantum depletion of the condensate. The main principles and a design 
of the detector of scattered atoms are suggested. 
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